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We have calculated the temperature dependence of the con- 
ductance variation (SS(T)) of mesoscopic superconductor nor- 
mal metal(S/N) structures, in the diffusive regime, analysing 
both weak and strong proximity effects. We show that in the 
case of a weak proximity effect there are two peaks in the de- 
pendence of 5S(T) on temperature. One of them (known from 
previous studies) corresponds to a temperature Ti of order of 
the Thouless energy (eTh), and another, newly predicted max- 
imum, occurs at a temperature T2 where the energy gap in 
the superconductor A(T2) is of order cth- In the limit L<f, < L 
the temperature T\ is determined by Dh/L^ (L^ is the phase 
breaking length), and not cth- We have also calculated the 
voltage dependence SS(V) for a S/F structure (F is a ferro- 
magnet) and predict non-monotonic behaviour at voltages of 
order the Zeeman splitting. 

Pacs numbers: 74.25.fy, 73.23.-b, 72.I0.-d, 72.I0Bg, 73.40Gk, 
74.50.+r 

Since the late 1970's it has been known that the con- 
ductance (G) of S/N mesoscopic structures depends on 
temperature (T) (and voltage (V)) in a non-monotonic 
way (see reviews [|]j|]). This behaviour was first pre- 
dicted in Ref. Q where a simple point S/N contact was 
analysed. The authors of Ref. ||, using a microscopic 
theory and assuming that the energy gap in the super- 
conductor (A) is much less than the Thouless energy 
ETh = fiD/L 2 (D is the diffusion constant), showed that 
the zero-bias conductance G coincides at zero tempera- 
ture with its normal state value (G„). With increasing 
T, G exhibits a non-monotonic behaviour, increasing to 
a maximum of G max « 1.25G„ at T m w A(T m ) and then 
decreasing to G„ for T > T m . 

Recently mesoscopic S/N structures have been fabri- 
cated in which the limit A >> exh is realised. In this 
case Nazarov and Stoof |J] (also see ) argued that the 
temperature dependence of the conductance G has a sim- 
ilar non-monotonic behaviour with a maximum at a tem- 
perature comparable with the Thouless energy, while si- 
multaniously Volkov, Allsopp and Lambert Q predicted 
that the voltage dependence of the conductance in an 
S/N mesoscopic structure (Andreev interferometer) has 
a similar form with a maximum at eV m ~ exh- This 
non-monotonic behaviour has been observed both in very 
short S/N contacts || and in longer mesoscopic S/N 
structures p0|- ^2| , pl| . In ref || it was noted that the 



conductance 5G = G — G n consists of two contributions. 
The first, SGdoSi is negative due to a proximity effect 
induced decrease in the density of states (DOS) of the 
normal wire which makes contact with a superconduct- 
ing strip |l3j . The other contribution SGmt (positive) is 
analogous to the Maki-Thompson (MT) contribution to 
the paraconductivity of S/N/S and N/S/N mesoscopic 
structures and was calculated in p^| . At T = and 
V = both contributions to the conductance are equal, 
as T or V increase the contribution SGmt dominates un- 
til a maximum is reached, then both these contributions 
decay. 

During the past decade a great deal of interest in 
the transport properties of N/S nanostructures has orig- 
inated from a desire to use superconductivity as a probe 
into phase-coherent transport. Indeed the above experi- 
ments reveal little about the superconductor itself, since 
the non-monotonic behaviour occurs at an energy much 
lower than A(0). In this paper we calculate the conduc- 
tance of mesoscopic S/N structures (see Fig. 1) over a 
wide temperature range (0 < T < T c ), and show that 
in the dependence of <5G(T), a second maximum may 
appear near T c when A(T m ) is of order the Thouless 
energy. Consequently, in contrast with non-monotonic 
phenomena studied to-date, this peak provides a novel 
quasi-particle transport probe for the energy gap of the 
superconductor. Indeed we show later that this second 
maximum is very sensitive to the damping rate inside the 
superconductor. We also show that if the depairing rate 
7 (for example, due to magnetic impurities) is not small 
compared to cth then the maxima in SG(T) occur at a 
temperature T\ « 7 and A(Ti) w 7. 

We consider S/N mesoscopic structures of the form 
shown in Fig. 1A and IB. Although they differ slightly 
from each other, in the limit li << 1 (l\ = L-y/L) the for- 
mulae for the conductances of these systems are identical. 
We assume the metals are diffusive and employ the well 
developed quasiclassical Green's function technique (see 
for example [lq] ) which has been widely used for study- 
ing transport phenomena in S/N mesoscopic structures 
1 iiH HflSil- Usin S the Keldysh formalism, the 
conductance variation of the structure shown in Fig. IB 
is given by @, 

5S = i°° d ^r%j w 
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where all the quantites are dimensionless (dimen- 
sional quantities will be denoted by a tilde); SS — 
(G~G n )/G ni [3 = (2T)-\Fl(e) = [cosh" 2 (e + eV) f3 
+ cosh~ 2 (e- eV)0\/2, t(e) =< tanh 2 (Re(u(e, x))) > 
, < A >= (1 — h)" 1 Jj dxA. All energies and voltages 
are measured in units of the Thouless energy (ct/i), the 
function u(e, x) is related to the condensate and nor- 
mal Green's functions: F R< ~ A ^ = sinh(w' R(j4) ), G R< ~ A ^ = 
cosh(u i? (' 4 )), which obey the Usadel equation, 



d 2 xx u R ^ - (k R ^f sinhu R ^ =0 



(2) 



where (k R ^) 2 = 7 =p lit. Eq. (Q) can be solved nu- 
merically JpLpl , and in some limiting cases analytically 

First we consider the simplest case, where the proxim- 
ity effect is weak, this occurs when the condensate func- 
tion in the N wire F R( - A ~> is small (\F R( - A 1\ w u R ( A ) « 1) 
HI). In this limit Eqs. ([!]) and (g) can be linearized. 
Also if the length L is longer than the phase breaking 
length (if 7 >> 1) then the N wire may be considered as 
infinitely long, thus the solution to Eq.(||) is determined 
by the expression p8[ , 



tanh(u/4) = tanh(uo/4)e" 



-hx 



(3) 



where uo is determined by the boundary condition at 
the S/N interface. In the case of a good contact at the 
N/N' interface the function u should be zero at this in- 
terface (x = 1), whereas Eq. (||) gives a nonzero value 
for u(e, 1). The correction to the solution (||) is small 
provided that 7 >> 1. In the case of the weak proxim- 
ity effect we solve the linearized form of Eq (Q) with the 
boundary condition at the S/N interface g], 

r b d x F R W = -F s R W\ x=0 (4) 

where r& is the ratio of the S/N interface resistance 
to the resistance of the N wire of length L; F^ = 
ia y F^ A ^ cos(</>/2) is the condensate Green's function in 
the superconductor, F* ' [A) = A/y/{e±iT) 2 - A 2 , T is 
the damping rate in the spectrum of the superconductor 
and <f> is the phase difference between the superconduc- 
tors. In the case of zero N/N' interface resistance the 
condensate function F R ( A > (or u^w) must equal zero 
when x — 1. The solution is 



pR(A) 



l(TyU 



R(A) 



F S R{A) sinh(fc fl ( A >(l - x)) 
rf,(fccoshfc) ii ( A ) 



(5) 



The proximity effect is weak provided that rbVT >> 1 
(a maximal value of F R ^ is achieved at e = A). In 
this case the function t(e) can be presented in the form 
%.t(e) = 1/4 < (F R - F A ) 2 >. Performing the spatial 
averaging we obtain in the limit l\ << 1, 



^) = (2^)- 2 ( Fco , 
V s r sin 
■ (fccoshfc) 2 L 2fe 



r sinh 2ki 
I 2fei 



sin 2k2 1 
2k 2 ' 



+Re( 



MM-l])cos 2 (</>/2) 
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where k = k\ + ife. The first term in Eq. (||) rep- 
resents the anomalous MT contribution to the conduc- 
tance variation and the second regular term is due to a 
DOS variation of the N wire caused by the proximity ef- 
fect. One can easily check that the partial conductance 
variation t(e) is zero at e = (in this case F R — F A ), 
increases as e 2 , reaching the first maximum when e is of 
order max[l,7] then decreases. At e m A the function 
t(e) has a second maximum (see solid line in Fig. 2). 
With the aid of Eq.(§) and @ we find the asymptotics 
of SS at low temperatures T << A (or (3 A >> 1). One 
has for the zero-bias SS at low temperatures, 



cos 2 (^/2) f Cl /3" 2 , 0» 1,(7 = 0) 
(2r 6 ) 2 C0 \5/7 7 / 2 /3 2 , 7- 1 «/?,(7»l) 



and at higher temperatures, 



(2r 6 ) 2 



c 2 , /3« 1,(7 = 0) 

c 37 - 1/2 , /3« 7 - 1 ,( 7 »l) 



(7) 



(8) 



where the coefficients are, cq rs 0.82, ci w 0.86, C2 ~ 
1.57 and c 3 w 1.56. It is clear from the expressions ([?] - 
|8|) for (J5 that the conductance variation SS has a first 
maximum at a temperature T\ rs max[eTh, 7] (we assume 
that both exh and 7 are smaller than the zero tempera- 
ture energy gap A(0)). 

Let us turn to the case of high temperatures T w A 
when (3 is close to the critical value (3 C = 3. 5/4 A (0). The 
contribution to SS caused by a variation in the DOS; cal- 
culated by summing over the Matsubara frequencies, is 
small and of the order SS « r^ 2 f3 3 / 2 (A(3) 2 . The main 
contribution is due to the MT term which can be pre- 
sented in the form, 



SSmt ~ Pc(2r b )-' 2 I de . 

lo ^/(e 2 - A 2 ) 2 + (2er) 2 



X(e) 



(9) 



Here and in what follows, we set <j> = for brevity. 
Where x(e) = |fccoshfc|- 2 [sinh(2A:i)/2fci -sin(2fc 2 )/2fc 2 ]- 
We replaced (3 in Eq. (|l|) by (3 C as j3 depends on A 
very weakly, (3 w (3 C [1 + 0.33(A/A(0)) 2 ] where the sec- 
ond term is very small (the calculations were carried out 
for A(0) = 20). For small energies (e << max [1,7]) 
the function x(e) is a constant equal to (2/3) when 
7 << 1 and equal to (7^ 3 ^ 2 ) when 7 >> 1, then 
for e >> max[l,7] x(e) decays as e~ 3 / 2 . The main 
contribution to SSmt stems from the singular region 
e ~ A, therefore in the main logarithmic approxima- 
tion we have SS M t ~ l/2^ c r 6 " 2 Aln(A/r)x(A) assuming 
that (A/r) >> 1. One can see that SSmt increases with 
increasing A from zero, reaches a maximum and then 
decreases. 

In Fig. 3 we present the dependence SS(T) calculated 
with the aid of Eqs (lil) and (0) . We see that besides the 
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main peak at T\ w 0.1T C (i.e. T\ w e^h), there is another, 
weaker peak in the conductance near T c . The tempera- 
ture T 2 at which the second peak is achieved corresponds 
to the condition A ss 7.6e-r/i- As the depairing rate 7 in- 
creases the first maximum is shifted towards higher tem- 
peratures, as predicted above. 

Fig. 4 shows the dependence SS(f3) for the same cross 
geometry in the case of a strong proximity effect (r& = 0), 
calculated using Eq. (||). We see that the weak peak 
near T c has disappeared and the position and height of 
the first peak depends on 7 essentially as before, moving 
to higher temperatures with increasing 7. 

In Fig. 5 the temperature dependence of the conduc- 
tance variation SS(@) is shown for the structure in Fig. 
1A. We assumed the weak proximity effect (r& = 5). The 
curves are presented for different resistances of the N/N' 
interface (r N / N i is the ratio of the N/N' interface resis- 
tance to the resistance of the N film). One can see that 
SS can be negative. The negative sign of SS is due to the 
shunting effect of the S strip which is stronger in the nor- 
mal state (in the superconducting state the S/N interface 
resistance increases) [p2| . As r N / N i increases, the height 
of the first maximum increases as the proximity effect is 
enhanced (the amplitude of the condensate function at 
the N/N' interface is not zero if rjy/jv 7^ and increases 
with increasing rjf/N')- The second peak in 5S(/3) be- 



comes weaker and disappears at 



N/N' 



= 10. This is in 



agreement with Fig. 4 where the strong proximity effect 
was considered (jb = 0). It is worth noting that in this 
geometry the second maximum is more pronounced (if 
r N/N' is not large) than in the geometry of Fig. IB .In 
Fig. 6 we show the dependence of the second peak on 
the damping rate T in the superconductor, as expected 
the second peak becomes less pronounced as T increases. 

In Fig. 2 we also plot the voltage dependence of the 
zero-temperature conductance variation SS(V) for the 
system shown in Fig. IB in which the N film is replaced 
by a ferromagnetic film (F). We assume a weak proxim- 
ity effect at very low temperature; in this case SS(V) « 
t(e) with e = eV. Fig 2 shows that the dependence 
8S(V) with increasing exchange h (measured in units 
tTh) changes drastically. First, the zero-temperature SSq 
is not zero at zero bias (time reversal symmetry is bro- 
ken) and has a non-monotonic behaviour as a function of 
h . Secondly, the low temperature peak is split and SS 
approaches zero atal/~/iif/i>> tTh- We note that 
these effects can be observed only in the case of a weak 
ferromagnetism when h (see for example Ref [^3| and 
references therein). 

In conclusion we have analysed the temperature de- 
pendence of the conductance S((3) for mesoscopic S/N 
structures of different geometries, and established that 
in the case of a weak proximity effect there are two 
peaks on the temperature dependence of 5S(T). The 
first maximum, as predicted earlier [^-|6|J^], corresponds 
to the temperature T\ of order of the Thouless energy, 



and the second corresponds to the temperature T 2 at 
which A(T 2 ) is of order tTh- Experimentally, this sec- 
ond maximum may already have been observed in fllpf (a) 
where a drop in the resistance R was observed near 
T c followed by a smooth increase in R with decreasing 
temperature (400mA" < T < 1200mA"). In contrast 
the first maximum occurs at a much lower temperature 
(7i 50 — 100mA). According to the theory presented 
above the position of the first maximum strongly depends 
on the depairing rate 7 in the normal wire, if 7 > tTh- 
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FIG. 4. Dependence of SS on the normalised inverse 
temperature, for the cross geomentry in the strong prox- 
imity limit. For different depairing rates 7. With 
A = 100,r = O.lJi = 0. 



(B) 



Superconductor 



FIG. 1. The structures considered. 
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FIG. 5. Dependence of SS on the normalised inverse tem- 
perature, for the geomentry shown in Fig. 1A in the weak 
proximity limit. For different interface resistances r NN i . 
With A = 10, 7 = 0.1,F = 0.1, h = 0.2. The inset is an 
enlargement of the second maximum around /3 » 0.1. Note 
SS has been scaled for convenience, by 0.1 for r NN i = 10,0.2 
for r NN i — 1 and 2 for r NN i — 0. 



FIG. 2. The dependence of i(e)(2r(,) 2 on energy, of the 
structure shown in Fig l.(B), for different values of the 
exchange field h — (solid line) , h = 1 (thick dot- 
ted line) and h — 7 (light dashed line), With 7 = 0.1, 
A = 20,r = 0.3,Zi = 0. 
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FIG. 3. Dependence of SS on temperature, for the cross ge- 
omentry in the weak proximity limit, for different depairing 
rates 7. With A = 20,r = 0.1, h = 0. The inset is an en- 
largement of the second maximum around T/T c — 1.0. Note 
the SS = 5S(2r b f. 



FIG. 6. Dependence of SS on temperature, for the geomen- 
try shown in Fig. IB in the weak proximity limit. For different 
damping rates F. With A = 20, 7 = 2,r = 0.1,h = 0. 
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